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Abstract. In order to stably grasp objects without using object models, tactile feedback from the 
fingers is sometimes necessary. This feedback can be used to adjust grasping forces to prevent a 
part from slipping from a hand. If the angle of force at the object finger contact can be determined, 
slip can be prevented by tire proper adjustment of finger forces. Another important tactile sensing 
task is finding tire edges and corners of an object, since they arc usually feasible grasping locations. 

This paper describes how this information can be extracted from the finger-object contact using 
strain sensors beneath a compliant skin. For determining contact forces, strain measurements are 
easier to use than the surface deformation profile. The finger is modelled as an infinite linear elastic 
half plane to predict the measured strain for several contact types and forces. The number of sen¬ 
sors required is less than has been proposed for other tactile recognition tasks. 

A rough upper bound on sensor density requirements for a specific depth is presented that is 
based on the frequency response of the elastic medium. The effects of different sensor stiffnesses on 
sensor performance are discussed. 
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1.0 Introduction 

Much of the recent work in tactile sensing has been devoted to recognizing objects and features. 
One method is to obtain an image-like array of an object profile using high density tactile sensors 
[Hillis, 1981; Overton and Williatnsm 1981]. This is useful for identifying the location, orientation 
and shape of an object with complicated surfaces, or identifying surface defects. Another approach 
to the recognition problem uses local low level tactile information and object models to recognize 
objects [Gaston and Lozano-Perez, 1983]. 

There are some low level tactile sensing operations that are useful for basic grasping, where the 
intent is keeping an object stably grasped in a hand, rather than recognizing it. The requirements 
for grasping polyhedra on a plane with two fingers without object models have been described else¬ 
where [Fearing, 1984], The most useful parameters to know are the surface normals, the angle and 
magnitude of a force at a contact, and whether the finger is touching a corner or edge. These 
parameters are a subset of those required to recognize general features. 

This paper attempts to show how this grasping information can be recovered using a simplified 
solid mechanics model and basic contact theory. Consideration of the mechanics has been infre¬ 
quent in the design and use of tactile sensors, (an exception is Kinoshita [1977]), but is important 
for relating sensory data to physical observables. There has been some work on sensing forces at 
fingers using tactile sensation to prevent slip [Stojiljkovic and Clot, 1977], Other techniques have 
relied on using rollers to detect the slip of an object in parallel jaw grippers [Masuda, Hasegawa, 
and Osako 1976], At a more advanced level, the tactile array approach will provide useful informa¬ 
tion for manipulation, such as finding a specific feature that is crucial to orienting a part accurately. 

A finger must have a compliant covering to take advantage of the increased prehension stability' 
possible at corners [Fearing, 1983], Another advantage of a soft skin is that contact areas are large 
enough to distinguish between features. For example, with a very hard skin, an edge and a side will 
contact with the finger at only a few points, and so may be indistinguishable without finger motion. 
To distinguish between them, it will be necessary to get either the displacement profile or tire con¬ 
tact stresses in a compliant finger covering. 

There are two different approaches to the tactile transduction problem. The first approach is 
based on measuring the deflection of a flexible membrane when it contacts a rigid object, or tire 
height of pegs touching an object [Page, Pugh, and Hcginbotham, 1976], The measurement is some¬ 
times done by an optical sensing scheme [Ozaki.Waku and Mohri, 1982], The second approach is 
based on forces beneath the surface changing electrical contact areas [Hillis, 1981], or compressing 
some resistive material without changing tire electrical contacts. 

We are interested in determining the actual contact stresses rather than just the contact profile. 
Although the profile aids in recognition of an object shape, it is only indirectly useful for determin¬ 
ing tire grasping forces. The problem is the indirect relation between the shape of the object and 
the profile of deformation, especially when there arc both normal and tangential forces applied at 
the contact. Because the compliant material is not perfectly compressible, it will tend to pile'up 
outside the contact region in a complicated manner. 

A mathematical problem with treating the contact as a deformation is that superposition does 
not hold. In general, forces superimpose, but displacements do not. For example, the stresses due 
to two points close together each indenting 1 mm. into an elastic medium are not equivalent to the 
sum of the stresses from each point indenting 1 mm. by itself. 

Consider a grain of sand pressed between a finger and a smooth surface. A profile change due to 
the height of the grain is very small, but high stress concentrations near the edges should be readily 
detectable. One researcher notes that a ridge of less than 20 pm height is perceptible for humans 
[Ledcrman, 1978]. 


2 



Fearing & Hollerbach 


Tactile Mechanics 


2.0 Using Strain for Tactile Sensing 

When determining contact stresses, there are practical reasons for measuring these stresses 
beneath the elastic material, rather than on the surface. Since we may be interested in both normal 
and tangential (friction) forces at the contact, we would need two different types of sensor at the 
surface. When sharp rigid objects indent an elastic material, very high stresses arc developed at the 
surface. These stresses are reduced by distance from the application area, so a fragile sensor would 
be better protected beneath a layer of skin. We note that biological tactile sensors are beneath the 
surface of the skin. A rather complete analytical and experimental study of stresses beneath the 
skin and tire mechanoreceptor responses to contacts without tangential forces was done by Phillips 
and Johnson [1981]. 

The general three dimensional analysis of stress and strain leads to very complicated expressions. 
If contact forces are constant along a finger, a two dimensional analysis will describe the behavior 
in a slice perpendicular to the finger axis. In the remainder of diis paper we analyze line, edge, and 
plane contacts for which a 2 dimensional analysis is exact for infinitely long contacts, such as a line 
load, and approximate for short contact lengths. We will speak of a point indenter or point force 
while meaning a planar section of a line contact. 

We now examine the 2 dimensional behavior of an infinite homogeneous elastic medium for the 
simplest contact case, a point force with negligible contact area. To find the stresses within the elas¬ 
tic half plane, we first find the stress distribution due to a concentrated normal force at the surface, 
as in Figure 1. The internal stresses due to a general contact stress on the surface can be found 
from the superposition principle. 

The analysis from Timoshenko [1951] shows that the internal stresses have a “simple radial dis¬ 
tribution”. That is, all the stress is in the direction of a radial line from the point of application of 
the force. 

From Timoshenko [1951], for the concentrated normal force we get 
— 2 P cos 9 


og = 0 


T r 6 = 0 (1) 

where a r is the radial stress at ( r , 6 ), og is the stress in the plane at ( r , 6) normal to the radial 
stress, T r g is the shearing stress in the r,8 plane, P is tire force per length, and r is the distance 
from the point of application. 

For a force inclined from the vertical by the angle a, 

_2 p 

o r = —- cos(a + 0) 

7rr 


og — 0 (2) 

Tr e = 0 
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Figure 1. Point Forces on Elastic Half Plane 


To get stresses in cartesian coordinates, we apply the tensor transformation: 
a x = a r cos 2 # 

Gy = a r sin 2 # 

T xy = o r sin#cos# (3) 

So finally, the internal stress due to a point force is: 
a x = —— cos(a + #)cos 2 # 

7rr 

— IP „ 

Gy = -— cos(a + #) sin 2 # 

* nrr 

_2 p 

t X y = -cos(a + #)sin#GOS# (4) 


where 


r = Jx 1 + y 1 

cos# = — , sin# = , cos(a + #) = —cos a - -^-sina 

r r r r 

All pressure sensors have an output based on die strain of the sensor, which is the fractional 
change in the linear dimensions of a small cubic volume clement. This volume change is to first 
order independent of the shear stresses, which change the angles of the faces in tire cube, and not 
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the dimensions. If we assume a linear elastic medium, we can apply Hooke’s Law: 


&X 



(5) 


where E is the modulus of elasticity (N/m 2 ), and e x is the strain along the x axis. 

For materials that are not completely compressible, a contraction along one dimension will be 
coupled with expansions along the other two. For an incompressible cube which undergoes a strain 
of 1% along one axis, there must be corresponding strains in the opposite sense of about .5% along 
the other two axes to maintain constant volume. The ratio of these two strains is called Poisson’s 
ratio, and characterizes the compressibility of the material. For a completely compressible material, 
Poisson’s ratio is 0. For a medium that does not change in volume with compression, such as a 
water filled sack or rubber, Poisson’s ratio can be taken as 0.5. Since rubber-like materials are 
popular for covering robot fingers and tactile sensors this assumption will be used. 

The relations between the stresses and strains are (Timoshenko and Goodier 1951): 


tx 

ty 


= Jy +ct z )| 

= ~jy jo^-Kerj+cO] 

= -jy [°t +°>)| 


( 6 ) 


where v is Poisson’s ratio = .5, and e is the strain along one dimension. 

There are two simplifying assumptions that can be applied to elasticity problems with the 
appropriate symmetry. For an infinite line force on an elastic half space, the strain in the direction 
of the line must be zero by symmetry. This is the plane strain assumption. Alternatively consider a 
slice of unit width out of the elastic half space, with a line load acting on it. On the faces of this 
slice, the stresses normal to the face must be zero to satisfy the boundary conditions. This is the 
plane stress assumption. 

For this analysis, we will assume a state of plane strain; there is no strain in the z direction. This 
assumption will be reasonable if the contacts are long compared to the finger radius. Phillips and 
Johnson [1981] found that an assumption of plane stress had a better qualitative agreement with the 
response of mechanoreceptors in finger skin, but the plane strain assumption will be used here 
because the contacts are assumed to be infinitely long. The form of the strains under these two 
assumptions is similar. 

So from the plane strain assumption e z — 0, 
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Figure 2. Strain for Line Load Normal to Surface 



Figure 3. Strain for Line Load Inclined 30 Degrees 
For the stresses due to the point indentor or line load (4), the strains from (7) become; 


for a = 0 


3 

-2 P 

4 E 

i nr 


cos 2 # — sin 2 #J 


cos(a + 6) 


3P | -22 

xcosa — ysina x —y 
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x ( x 2 -y 2 ) 


( 8 ) 
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Figure 4. Strain Sensor Based on Resistance 


Figures 2 and 3 show the strains in the homogeneous elastic medium for a line load normal to the 
surface and one inclined 30 degrees. 

2.1 Measuring Strain 

Since our plane strain assumption gives a change in volume that depends only on the strain in 
the x direction, there are some implications for volume dependent strain sensors. Figure 4 shows 
how a transducer based on the volume resistance could be used to determine the strain at a certain 
depth. Consider a block of conductive material with the electrodes on the top and bottom. The 
resistance is given by: 

R _ £k _ ££ ( 9 ) 

A yz 

where p is the resistivity, R is the resistance, L is the length, and A is the area. 

For plane strain the new dimensions are: 

x' = Jf(l + £ x ) 

y' — ^(1 + e^,) = y(l-e x ) since e y = —s x 

t 

z — z 

For a unit volume clement 

R = '^_ e ) ~ p(l + e .v)U+ £* + £.?+ • • • ) ~ p(l +2e v ) (10) 


This transducer would output the variable of interest. A similar result can be shown for the plane 
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Figure 5. Three Equally Spaced Strain Sensors 


stress assumption: 

R = (l + e y Xl + c,) ~ P( l + e x~ £ y~ e z) = Pd + 2s x ) 

since -e z = e^ + e, (11) 

3.0 Determine Force Angle for a Line Load 


Can the angle of inclination, location, and magnitude of a line load be recovered from die type 
of sensor described previously, that measures only one parameter, e v ? Assuming negligible sensor 
dimensions, it is straight forward to set up the equations for three strain sensors equally spaced on a 
horizontal plane beneath the surface. Referring to Figure 5 they are: 



3 P |.xcosa — (y + Zf)sina 
2 t jE | x 2 + {y + b) 2 


iy + b)‘ 


— 3Ppc cos«—ysina 
2 t TE(x 2 +y 2 ) 2 " 



(12) 


-3P Lvcosa- 


-(y — Zt)sin«J 


2tjE\x 2 + (y — b)‘ 


iy — b)‘ 


where x is the sensor depth and b is the sensor spacing. 
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««, = -18 
= 19 


Figure 6. Finding the Force Location with Multiple Roots 


We can use “equation counting” methods [Richards, Rubin, and Hoffman, 1981] to determine if 
these 3 measurements arc sufficient to uniquely determine the line load on the surface. The impli¬ 
cit function theorem [Sokolnikoff and Redhcffer 1966] states that if tire system of equations (12) is 
continuous and has continuous first partial derivatives with respect to the independent variables, 
and if the Jacobian is non-zero, then there will be a unique and continuous solution to those equa¬ 
tions. The determinant of the Jacobian of eq. (12) is not identically zero, but does disappear for 
some values of parameters. 

To determine where the extra solution points are, tire system of equations was solved by elimi¬ 
nation of the magnitude and angle of the force, to get a polynomial in y and the measured strains: 

/O') = 2e Xj [x 2 + y 2 } [x 2 - O'-*) 2 ] [x 2 -(y +Z>) 2 ] 

- e Xl [x 2 + (y + 6) 2 ] [x 2 - y 2 J [jc 2 - 0 — A-) 2 ] (13) 

- £ X3 [x 2 + (y-fc) 2 ] [x 2 ~>' 2 ] [x 2 - (y + Z>) 2 ] = 0 
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SIDE VERTEX 



FINGER 


a) Contacts with only normal force 



c) Contacts with generalized forces 


Figure 7. Finger Contacts with Side and Vertex 


This eighth degree polynomial has multiple real roots close to the sensor locations. It is interest¬ 
ing to note that if the sensor depth and the spacing are equal (x = b), one of the false roots will 
always be at y = 0. Constraints on the possible forces, such as magnitude ranges and angle limits, 
allow some of these solutions to be rejected. For example, a reasonable constraint is that all forces 
are directed into the surface (compressive). Tension will not occur at the contact without an 
adhesive bond. If there are still possible roots left, a fourth sensor measurement will be necessary 
to get a unique solution. 

Figure 6 shows an example with false roots that can be discarded. Two can be rejected because 
the point force would have to be in tension, and one other because the force is too close to hor¬ 
izontal and would slip. 


4.0 Distinguish a Vertex from a Side 

In this section, we will try to differentiate between the strains beneath the skin due to the vertex 
and tire side contacts. Figure 7 shows the assumptions used for these contacts. Assume the side is 
perpendicular to the finger and is moving perpendicularly into the skin. Foi tire vertex, assume its 
center line is perpendicular to the edge of the finger and is moving perpendicularly into the skin. 
These can be viewed as contacts with a normal force but no moment (Figure 7a). 
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Figure 8. Two Extreme Vertex Cases 


Figure 7b shows normal contacts with a tangential force. The most general case of a moment 
and a traction force is shown in Figure 7c. Restricting the analysis to the 2 dimensional case with 
plane strain, wc will attempt to characterize the strains that would be measured for normal contacts 
without a tangential force. 

For this analysis, we will ignore traction from friction while indenting normally into the elastic 
medium. With friction present, the elastic medium would adhere to the object and cause surface 
shear stresses as it tried to flow away from the indentation. The normal stresses are not affected 
very much by this traction effect [Engel and Conway, 1971]. 

What does a vertex feel like? Consider die vertexes touching die finger as in Figure 8. The first 
case appears to be the right side of a rectangular indenter. If we could consider the finger and the 
contacting side extending far away from the right side, then the left side of the indenter would have 
no effect on the right side. The surface stress for a rectangular indenter on an elastic half plane is 
given by [Conway, Vogel, Famham and So, 1966]: 


CTjCx = 0) = 



for \y | <a 
otherwise 


(14) 


where P is the force per unit length, and a is the half width of the contact. 

The second case is equivalent to a wedge of 90 degrees being pressed into an elastic half plane. 
This stress is given by [Gladwell 1980: Robinson and Thompson 1974] as: 


a x {x = 0) = 


Kcol—cosh 1 -~ 
■2 IT 

0 


for | y |<a 
otherwise 


(15) 


where K is a function of Poisson’s ratio and the modulus of elasticity, and 6 is die vertex angle. 
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Block Indenter 



Wedge Indenter 



Figure 9. Stresses and Strains for Vertexes 


Figure 9 shows the stresses and strains for these vertex contacts. A vertex looks very much like a 
line load especially for small indentations and small vertex angles. If the line load and vertex have 
approximately the same strain profiles, the method of section 3.0 could be used to find tire angle of 
force at this contact. 

Now consider the contact of a stiff flat side and a circular elastic finger as in Figure 10a. This 
two dimensional problem is difficult to analyze analytically, so we will apply Saint-Venant’s Princi¬ 
ple, as explained by Timoshenko [1951]: 

This principle states that if the forces acting on a small portion of the surface of an elastic body are 
replaced by another statically equivalent system of forces acting on the same portion of the surface, 
this redistribution of loading produces substantial changes in the stresses locally but has a negligible 
effect on the stresses at distances which are large in comparison with the linear dimensions of the sur¬ 
face on which the forces are changed. 

If the area of contact is small compared to the finger radius, we can approximate the area near 
the contact as an infinite elastic half plane. This approximation will be most accurate near the 
outer surface of the finger in the quadrant where the contact is. It is also necessary to reformulate 
the contact stresses for the cylindrical finger to match the half plane approximation. 

In the solid mechanics literature, the common problem is the rigid cylinder indenting an elastic 
layer, the punch problem, rather than an elastic cylinder contacting a rigid plane. Bental'l and 
Johnson [1968] have done an analysis where the elasticity of the cylinder and layer are equal. 
Hahn and Levinson [1974] have analyzed the stresses for the contact of a rigid cylinder and a 
cylinder covered with an clastic layer. That analysis is probably most relevant to the rigid finger 
with an elastic covering contacting a rigid surface, but as in much of the contact literature, the form 
of die solutions (infinite series) is inconvenient to work with. 

For a rigid cylinder indenting an elastic half plane, as in Figure 10b, the stress on die surface is 
given by Conway [1966] as: 

o x iy) = fa 1 -/ (16) 

7 red v 

where a is the half widdi of the contact region and P is die force per unit length on the cylinder. 
The stress and strain for this contact arc shown in Figure 11. 
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-a +a elastic half plane 


b) 


Figure 10. Approximating a Contact with a Side 



-10 -5 0 5 10 

Figure 11. Stresses and Strains for Cylindrical Indentor 


The half width of the contact from [Bentall and Johnson, 1968] is: 

l 


where R is the radius of the cylinder. 

For equilibrium in the contact region, the stresses should be equal and of opposite sign on the 
cylinder and the surface. Now for an elastic cylindrical finger contacting a side, there should be no 
infinite stresses such as occur at shaip edges of indentations. Also, the stress will have a peak value 
in the center of the contact, and approach zero at the edges. Using this justification, we will assume 
that the stress on the finger can be approximated by cq. (16) for small indentations. 


a = 2 


4RP(l — v 2 ) 
•uE 
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Contact Stress 



Figure 12. Tangential Stress Proportional to Normal Stress 


The stresses due to contacting a side and a vertex are distinguishably different. For the side, 
die re is a smooth, continuous stress function. For the vertex, there is a discontinuity in the stress at 
the tip of the wedge. (Actually, the medium deforms around diat discontinuity, spreading the stress 
over a finite area). It should be easy to distinguish the two just by looking at the width of die strain 
pattern. For the same elasticity and pressure, the wedge has greater indentation depth with a 
smaller contact area than the side contact. If the vertex angle is less than 90 degrees, there will be 
a significant difference between the two contact widths. Arguably, diis task would be simpler if the 
indentation profile were available. 

5.0 Determining the Side-Finger Contact Force 


For a finger touching the side, we will assume that die contact is of the form shown in 
Figure 7b, widi bodi normal and tangential components, but no moment. The tangential com¬ 
ponents will be modeled as all being in die same direction widi a force that is directly proportional 
to the magnitude of die normal force [Smith and Liu, 1953]. Figure 12 shows this assumption for a 
general pressure pattern. (The tangential force must be less than pN, the friction force). 

A very complicated expression for the stress below the surface for an elliptical stress distribu¬ 
tion on the boundary was found by Smith and Liu [1953]. Wc do not need the complete stress 
expression, but want to determine the angle of force at the contact. For grasping applications, the 
angle of force a is useful for predicting slip. 


To find the strain underneath for the stress for the cylindrical indentor from eq. (16), we 
integrate along the surface where 


P(y) = 


IP 

2 / a l + a 2^ 

a \ (y 0 ) 

7ra z 

Z 

[ 0 otherwise 


for a x < y < a 2 


(18) 


where a x and a 2 are the edges of the contact, and a = a 2 —a x . Applying the superposition 
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Relative Strain at Depth 
2.0 for Contact Width 10.0 




oC= 0.3 
oC= 0.6 
oC= 0.9 


Figure 13. Strain for Cylinder Contact with Tangential toads 


principle to the strain due to a line contact cq. (8): 


oo 


e*U^)= / 
-00 


-3 jrcosa - j- o sin«J|r 2 -y’ o 2 j/ > 0' -y 0 ) 


2-nEr A 


dy 0 


(19) 


Substituting the integration limits for P(y) gives: 


e x (x,y) = 


-3 P 
•n 2 a 2 E 


y~ a \ 


/ 

y-a 2 


x cos« — y 0 since 
(x 2 +y 0 2 ) 2 


-y 0 


■(y-y„- 


a x + a 2 


’ dy 0 


( 20 ) 


where r 2 = x 2 +y 0 2 . 

Figure 13 shows the strain for the cylindrical indentor with various tangential loads applied. 
Note that the slope of the strain underneath the center of the contact increases with increasing 
tangential force. The slope of e x (x,y) at y = (a\ + ai)/2 can be found by evaluating the derivative 
of eq. (20), using the first fundamental theorem of integral calculus to get: 


de x 
dy 


3/3 P 
2nr 2 aE 



„ 2 

r^l 

* ,+ «f > 


sina 


(21) 


When the contact length is large compared to the sensor depth, this simplifies to 
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Surface Stress Elastic Medium 


Discrete Strain Measurements 



e x (nb) 


Sensor Spacing b 


Figure 14. Discrete Strain Measurements 


de x _ — 6\/3/ } sina ,~ 2 ) 

dy ~ 7 t 2 o 3 E K ’ 

So for a given contact, a can be easily recovered if the pressure can be estimated. 

The more complicated contact with generalized forces of Figure 7c can be treated as the super¬ 
position of stresses due to forces plus the stresses due to the moment load. The elasticity literature 
deals with moments applied to punches indenting an elastic medium [Alblas and Kuipers, 1970; 
Galin 1961]. In those cases, the pressure distribution due to die moment load is an odd function 
symmetric about the origin, with half the surface in tension and the other half in compression, so 
the moment can significantly affect the internal strains. 

6.0 Sensor Density Requirements 

In section 3, a method to determine the locadon, direction, and magnitude of a line contact was 
oudined that required at least three strain sensors. A more general question is, how many strain 
sensors are needed to identify a general stress pattern on the surface? This problem is considered 
in the context of the spatial frequencies of the applied stress. The spatial frequency of the stress 
pattern is just die variations in pressure per unit length, not time dependent variations. We shall 
assume that die pressure pattern is approximately bandlimited, that is, most of die stress spectrum is 
concentrated at low frequencies. This is the case for contact stresses in the form, of (15) and (16). 

Figure 14 shows a linear system representation of this discrete measurement process, where the 
sensors are of negligible width. The sampling theorem says that any bandlimited signal can be 
recovered from its discrete samples if the samples are taken at a sufficiendy high frequency, the 
Nyquist rate. (For example see [Oppenhcim and Schafer, 1975]). Here the sampling frequency is 
the minimum strain sensor spacing needed to avoid aliasing when recovering the continuous strain 
measurements from the discrete samples. If there are not enough sensors, die strain due to the high 
frequency components will appear as noise when the continuous strain is recovered from the sam¬ 
ples by low pass filtering. 
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Figure 15. Spatial Frequency Response of Elastic Medium 



Since the strain at a point beneath the surface is a linear function of the two stress components 
a x and a y for plane strain, the total strain at a point in the elastic medium can be found by the 
superposition of the strains underneath due to each stress component on the surface. It should be 
possible to treat the normal and tangential stresses separately, and determine the highest frequencies 
of interest for each. If the force on the boundary has no tangential components, the pressure distri¬ 
bution Piy) will be a scalar at each point. The supeiposition integral is: 

00 

t x {d,y 0 ) = J ~ip-h(yo-y) dy (23) 

-00 b 


where Piy) is the pressure distribution on the surface and h(y) is the strain at depth d due to a 
unit normal pressure point at the origin. From equation (8): 


h(y) = 


-3 d(d 2 ~y 2 ) 
2tt {d 2 +y 2 f 


(24) 


The spectrum of the “impulse response” h O') can be found from simple properties of Fourier 
transforms: 


t j(d 2 +y 2 ) 
d 2 

t r(d 2 +y 2 ) 2 


e -d\2m\ 


7re 


-d 2irs 


din 


+ s 


for s > 0 


df. 


v (d 2 +y 2 ) 2 - 


•nde 


-dlTTS 


dim 


for s > 0 
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Strain Without Ridges Strain With Ridges 



Figure 16. Effect of Ridges on Strain Amplitude 


— 3 d(d 2 —y 2 ) 
lAcP+y 1 ) 2 


— 3 •ndse~ d \ 2 ' ns \ = H(s) 


(25) 


The incompressible medium in plane strain has no response to a zero spatial frequency, that is, 
a constant pressure over an infinite extent, because there is no space for the medium to escape into.’ 

The frequency response of the elastic medium as shown in Figure 15 is not strictly bandlimited, 
but it does have very steep skirts. We shall assume that the overall frequency response P(s)H(s ) 
has a high frequency response like H(s), where P(s) is the pressure spectrum. 

We can choose a sampling frequency by deciding how much aliasing is permissible. This 
requires an engineering judgement based on the estimated accuracy of the strain sensors. It seems 
that any aliasing that was less than one part in 1000 of full scale would be negligible because most 
sensors would not have a dynamic range much better than that. Also, non-linearities and tempera¬ 
ture instabilities get to be of at least the same order of magnitude. 

The following rough estimate can be used for sampling requirements. We want the value at the 
tail to be down a factor of 1000 from the maximum which occurs at lirsd = 1. So: 

-~y = 10 3 , or / = 10.2 (26) 

The samples should be at twice f so 2rrsd > 20.4, or sd > 3.2. Here s is tire number of samples 
per unit length. 
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For a sensor depth of 1.0, a sensor density of 3.2 sensors per unit distance would allow recovery 
of the continuous strain response due to a normally directed surface stress distribution with aliasing 
errors down at least a factor of 1000 from the desired signal. This sampling rate is probably too 
conservative. The types of features touched may have rounded edges, and start out with a more 
band limited signal, so a lower sensor density could be adequate. The effective sensor density can 
be increased by multiple measurements as a finger is moved along a surface, so the uldmate fre¬ 
quency response is limited by the number of samples per second, and the sensor accuracy. 

6.1 Using Finger Ridges to Enhance Sensitivity 

The infinite elastic medium has a peak response when the wavelength equals the strain sensor 
depth times 2tt. There seems to be an interesting possibility of maximizing die strain sensors’ out¬ 
puts by locating them at the maxima and minima of strain, and predetermining the frequency and 
phase of the pressure. This could be done in principle by inserting a thin and flexible sinusoidal 
grating between the finger and the object touched, which could superimpose a sinusoidal stress on 
the regular stress. A simple method of adding the sinusoid of the desired frequency is to use a 
ridged finger. 

In Figure 16 we compare the strain amplitude between the unridged and ridged fingers, where 
the ridges are of the optimal wavelength. In both cases the sensors are at the same depth. In the 
unridged case the total contact stress is about 1.5 times the ridged case. Even neglecting the 
difference in total contact stress, the ridges can enhance the amplitude by a factor of 2. The sensors 
should be located beneath and between ridges to detect the maximum amplitude peaks. This tech¬ 
nique is similar to the electrical chopper that allows a DC level to pass through an AC coupled cir¬ 
cuit. 

It is curious to note that in the human finger, there is a mcchanoreceptor located directly 
beneath every papillary ridge [Lederman, 1978]. Phillips and Johnson [1981] gave the depth of this 
strain measuring mechanoreceptor in the macaque monkey as about 500 pm. This is about the 
same depth as in the human finger [Quilliam 1978]. The optimum ridge spacing for this depth is 
about 3 mm, much larger than the papillary ridge spacing on the authors’ fingers. However, it is 
interesting to speculate that fingerprints may still have some role in enhancing strain measurements. 


7.0 Discussion 

What effect does the skin compliance have on tactile sensing abilities? A stiffer skin will 
develop greater contact stresses for the same indentation, but the subsurface deformations (the 
measured strain) will be the same [Phillips and Johnson, 1981]. The smaller contact area and defor¬ 
mation with a less compliant material helps to make the approximations we have used more valid. 
Greater stiffness in the skin may also serve to protect the sensors because a stiffer skin has more 
resistance to puncture. 

A very compliant skin on a sensor may have saturation effects sooner than a harder skin. Fig¬ 
ure 17 shows three fingers applied to a surface with the same normal force, but with different cov¬ 
ering softness. For the very stiff skin, the contact is not large enough to determine the size of the 
particle. For the very soft skin, the stress due to the small particle may become proportionally less 
significant titan the stress due to contact with the surface around the particle. One way to choose a 
skin compliance for good sensing characteristics would be to decide on a typical finger force range 
and the smallest contact area of interest, and use this information to determine the necessary 
modulus of elasticity. 

An area for further research is tire spatial bandwidth that is required for feature recognition. 
The two-point discrimination test [Valbo and Johanson 1978], seems to be tire classic test for tactile 
sensor resolution. This test would seem to require a very high bandwidth because of the two 
impulse functions, but it may turn out that only lire low frequencies are important. While we are 
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Figure 17. Contacts for Different Skin Compliances 


not interested in resolving between two point indenters. this problem may be useful for setting a 
smaller upper bound on the density of sensors. 

In psychophysical experiments, subjects arc instructed that there will cither be one or two points 
indenting their skin. It appears that the sensors should be able to distinguish between these two 
eases for two-point distances significantly smaller than the individual receptive fields of each sensor. 
Here, the resolution is not limited by the density of the sensors, but by the measurement accuracy 
of the sensors. 

it is useful to consider each strain sensor as having a response function that overlaps its neigh¬ 
bor. A point force on the surface will cause an output for each sensor that is proportional to die 
height of the response function underneath that point. For two point forces on the surface, the out¬ 
puts will be the linear superposition of the responses of each point force. Wc can measure the out¬ 
put in each sensor channel, and determine if it is consistent with just a single point force. How¬ 
ever, there might be different combinations of point forces that can give the same output. A very 
loose analogy can be made to the color matching problem [Richards, 1979], where several different 
stimulation combinations can give a white sensation. 

In Figure 18 there are locations where two point forces can be applied that will give a response 
in the two sensors equivalent to a single point force. In Figure 19, a third sensor has been added to 
allow discrimination between all the two point and one point cases, limited by the measurement 
accuracy. It will be difficult to distinguish points that give responses far down on the “tails” of the 
channel. 

It is probably reasonable not to spread the sensors too far apart. A guess is that the practical 
maximum separation would be less than twice the depth of the sensors. (Note that for the two- 
point discrimination task, the modulus of elasticity has no effect on the width of the strain response, 
so it will have no effect on the resolution). More work should be done on choosing the depth and 
sensor spacing to optimize discrimination. 

It seems that some recognition tasks would be easier using a sensor that responded to deforma¬ 
tion. That would make distinguishing a vertex and a side easy work. However, other tasks, such as 
sensing imminent slip, appear easier with the strain sensor approach. Sensors that combine surface 
shear sensors, depth strain sensor’s, and surface deflection sensing would simplify the problem con¬ 
siderably. 


20 




Fearing & Iloilerbach 


Tactile Mechanics 



Figure 18. Two Strain Response Sensors 



Figure 19. Three Strain Response Sensors 


8.0 Conclusion 

In this paper, we used a simple linear elastic model to predict strains beneath a compliant skin 
for a finger touching a knife edge (line load), a corner, and a flat surface. Appropriate approxima¬ 
tions were made to get simple expressions for tire contact stresses for a cylindrical finger contacting 
a side. In general, three strain measurements were shown to be necessary for determining the loca¬ 
tion, magnitude, and direction of force for a line contact, but degenerate cases with multiple solu¬ 
tions arc possible. The angle of force was shown to be related to tire slope of the strain curve for a 
contact with a side. Further work needs to be done in recovering the angle of force from contacts 
with other curvatures. 

The elastic medium was examined as a signal processing element, which led to tine interesting 
possibility of maximizing strain sensor output by adding a grating to tire sensor surface with a 
period equal to the sensor depth times 2tt. 
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We have not attempted to determine stresses and strains for the generalized three dimensional 
contact. At least five strain measurements arc required to detennine the magnitude, two angles, and 
y,z location of a point force on an clastic half-space. These equations are quite complicated [Mind- 
lin, 1936]. It may be possible to find this information by combining two orthogonal planar solu¬ 
tions, but this was not attempted here. 

The discussion here seems to indicate that four sensors, with fairly high dynamic range, may be 
adequate for the planar case for finding force magnitudes and directions, and for distinguishing 
between two contact types (vertex and side). With a thick skin, one can get a lot from a few sen¬ 
sors. This number of sensors will probably not suffice for determining curvature and more compli¬ 
cated force resolving problems. 
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